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Abstract 

Recently, a new technique for calculating the defect energetics of alloys based on Equiv- 
alent Crystal Theory has been developed. This new technique successfully predicts the 
bulk properties for binary alloys as well as segregation energies in the dilute limit. We 
apply this method for the calculation of energy and force as a function of separation of 
an AFM tip and substrate. This study was done for different combinations of tip and 
sample materials. The validity of the universality discovered for same metal interfaces 
is examined for the case of different metal interactions. 
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The atomic force microscope (AFM) [1] has enabled the measurement of interfacial 
forces at the atomic level as a function of separation [2], However, the ability to perform 
theoretical calculations of these forces is limited by the complexity of tip geometries. Also, 
the need to calculate force as a function of separation and the need to include relaxation 
makes first-principles calculations difficult to perform. A new formalism, Equivalent Crystal 
Theory (ECT) [3], which has given accurate surface energies and relaxations for fee, bcc 
metals and semiconductors [4] has been used to calculate the force versus separation for a 
range of tip geometries [5]. However, ECT has only been formulated for single component 
materials, whereas it is necessary to be able to calculate the same properties for different 
materials in contact. 

Recently, a new method for treating multicomponent systems, which builds on the ideas 
of ECT, was introduced by Bozzolo, Ferrante and Smith (BFS) [6]. Here, we give a brief 
description of the method and present an application to a number of fee and bcc metals in 
contact for several tip geometries. We will present the calculations of the binding energy 
and force as a function of separation and discuss the universal aspects of these curves. 

In BFS, the total energy of a multicomponent system is the sum of the energies of 
individual atoms. For each atom, we break up the energy into two parts: a strain energy 
e s and a chemical contribution e^, linked by a coupling factor g: 

e * = e ? + 9i{e? — € o,) ( 1 ) 

where i denotes the atomic species of a given atom. 

The strain energy, ef , is obtained from a standard ECT defect calculation where all 
the atoms surrounding the atom in question are taken to be of the same atomic species i. 
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In the spirit of ECT, an equivalent lattice parameter of is straightforwardly computed and 
the strain contribution becomes 

ef = E' c r(af *) ( 2 ) 

where F*(a) = 1 - (1 + a)e~ a and E' c is the cohesive energy of an atom of species i. The 
argument of* is given by of* = (of - o‘)/fc, where o' and Z, are the lattice parameter and 
screening length of a pure crystal of species i. The coupling term gi is defined in terms of 

the scaled lattice parameter of* a s <7; = exp(—af*). 

The chemical contribution ef is obtained by yet another standard ECT calculation. 
As opposed to the strain energy term, the surrounding atoms retain their chemical iden- 
tity, but are forced to be in equilibrium lattice sites. If Ni k denotes the number of 

nearest(next)-neighbors of species k of the atom in question (of species i) then the ECT 
equation [4] to be solved for the equivalent lattice parameter is 

NR\'e~ 0iR ' + MR%e- (a '+*l )R2 = £ N ik r\'e~ a " T ' + £ M * r 2 ' r ( °“ + £ )rj ( 3 ) 

k k 

where N(M) is the number of nearest (next)- neighbors in the equivalent crystal of species 
t and Ri(R 2 ) is the nearest(next)-neighbor distance in the equivalent crystal of lattice 
parameter af . rj and r 2 , are the equilibrium nearest- and next-nearest-neighbor distances 
in an equilibrium crystal of species t, respectively. The chemical energy is then computed 
with 

ef = 7 £^*(of*) (4) 

and 

ef = 70 E' c F'(aC;) (5) 
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where 7(70) = +1 if af*(a£*) > 0 and 7(70) = -1 otherwise, and af‘ = (af - a*)/£. The 
scaled lattice parameter af* is obtained from Eq.(3) with the parameters a,* listed in Ref. 
6, and af * is computed by solving Eq.(3) but with a,-* = a,-. The rest of the parameters 
appearing in Eq.(3) are listed in Ref. 6. 

An overview of typical results is shown in Figs. 1-3. Fig. 1 displays the comparison 
between the experimental values of the heat of formation [7] of selected binary compounds 
A x Bi- x and the predictions of this method, for the full range of concentration x. In Fig. 2 
we compare the predictions for the lattice parameters of the binary alloys of Ni, Cu and Pd 
with experimental values [8]. The method’s sensitivity is highlighted in Fig. 3, where the 
lattice parameters of Ag-Au alloys are compared with experiment [8]. A lattice contraction 
is experimentally detected for certain concentrations involving significantly small changes 
of the lattice spacing, and they are successfully reproduced by our method. Similar results 
are obtained for bcc alloys. A full description of the method is given in Ref. 6. Even though 
BFS is a semiempirical method, its dependence on experimental input is minimal in that 
only two experimental values (i.e., the heats of solution in the dilute limit) are needed and 
most of the formalism is based on the behavior of the pure elements. 

The only limitation is the availability of reliable experimental data for the correspond- 
ing alloys. BFS is designed to provide, with minimum computational effort, a systematic 
description of arbitrary systems of more than one atomic species. Because of its numerical 
simplicity, this method, together with a Monte Carlo procedure to account for individ- 
ual atomic displacements, is particularly well suited for computer simulations of typical 
STM/AFM experiments. 

In this first application of BFS to the study of metallic adhesion using AFM, we consid- 


4 



ered several bimetallic junctions, concentrating on the effect of the tip shape on the adhesive 
forces in AFM. We left out the possibility of individual atomic displacements and avalanche 
effects for small gap widths, as our concern is with the existence of a universal behavior 
that include bimetallic systems in ideal (i.e., rigid) conditions. If these effects were taken 
into account, one would expect a greater sensitivity of the energetics of the tip geometry as 
well as tip and sample materials. We restricted our calculations to systems in which both 
the sample and the tip are fee or bcc metals. Although there is some experimental evidence 
that multiple atom tips are typical in the scanning process [10], in this initial study we 
choose to model the interaction between the tip and the substrate with a single-atom tip. 
We will show that different tip geometries account for minor corrections to the energetics 
of the system. This result is consistent with the body of available STM/AFM data [10]. 

Fig. 4. a shows the variation of binding energy with separation for two different 
STM/AFM A1 tips on a (100) A1 sample. The tip-substrate separation is taken to be 
zero when the tip atom is in contact with the sample surface (in this calculation, this corre- 
sponds to the fcc-hollow site). The two tips examined are sketched in Figs. 4.b and 4.c. As 
discussed before, the differences in energy due to different tip geometries are rather small. 
Moreover, the adhesive energy (i.e., the minimum of the binding energy curve), if amenable 
to measurement, should give an indication of the tip shape. 

Fig. 5 displays the dependence of the binding energy curve in terms of lattice symme- 
tries. Two different surfaces of Ir are considered: (100) and (111)- As expected, the adhesive 
energy values are smaller for the (111) face. These curves correspond to a pyramidal tip 
(see Fig. 4.c), but the essential features should be the same for other tip shapes. 

Fig. 6. a and 6.b summarize the binding energy versus separation results for a variety 
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of fee and bee materials (Al, Ni, Cr, Ir) for their (100) face and a pyramidal tip. Consistent 
with previous findings, all curves scale onto one curve which is accurately described by a 
Rydberg function. Similar conclusions can be drawn for the force vs. separation curves, 
shown in Fig. 6.c. 

Finally, Fig. 7 displays the variation of binding energy with separation for a variety of 
tip-substrate systems A-B (A:tip-B:substrate) where Ni is one of the components. There is 
a noticeable change in the qualitative behavior of these results as the lattice parameter mis- 
match between the components changes. However, we found that for bimetallic junctions, 
the binding energy curves also scale onto a unique curve which can be fitted satisfacto- 
rily Rydberg function,’ as in the case of pure metals as shown in Fig. 7.d. This result is 
independent of lattice structure, crystal face, tip geometry and materials involved. 

Simulations like the ones presented in this work are computationally extremely simple, 
allowing for straightforward extensions of these ideal situations, including multiatomic tip 
effects and, more importantly, tip induced relaxations. In particular, BFS, together with 
the recently developed theory of avalanche in adhesion [11], should be ideally suited for the 
study of stability of a tip and sample at very small separations [12]. 

Conclusions 

The semiempirical method for the treatment of multicomponent systems developed by 
Bozzolo, Ferrante and Smith is physically accurate and computationally simple. Because of 
the ability of this method to effectively deal with systems with more than one atomic species, 
it provides an ideal framework for theoretical modelling of defect structures such as impurity 
atoms. For the case of STM/AFM, applications to a variety of tip-substrate combinations, 
together with some insight on the dependence of the tip structure yield encouraging results, 
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in the sense that in all the situations considered, the universal behavior expected (and 
observed) for pure metals is also reproduced. 
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Figure 2. — Comparison between the lattice parameters obtained with BFS (solid curve) and Vegarti's 
law results (dashed lines) for different fee binary alloys. Experimental values are shown for Nl-Cu 
(solid squares), Nf-Pd (solid circles) and Pd-Cu (solid triangles). 
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TIP-SAMPLE SEPARATION [A] 

(a) Variation of binding energy wtth separation for a STM/AFM Al tip at the fcc-hotJow site of a At (1 00) sur- 
face for a one-atom tip (solid curve) and a pyramidal tip (dashed curve). The separation is taken to be 


zero when tip and \ 
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Figure 4. 
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Figure 5.— Variation of binding energy with separation for an Ir tip on a (100) surface (solid curve) 
and a (1 1 1 ) (clashed curve) Ir surface. The separation is taken to be zero when tip and substrate 
are In contact 
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TIP -SAMPLE SEPARATION [A] TIP-SAMPLE SEPARATION [A] 

(a) Variation of binding energy with separation at the hollow site on (b) Scaled binding energies end separations fot the cases shown 
a (100) surface of different fee and bcc metals: Al/AI (short dash), In (a). 

Nl/Nl (dots), Cr/Cr (long dash), Ir/lr (solid curve). 



tip-sample separation [a] 


(c) Variation of the force with separation for the cases shown In (a). 
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(c) Nl/Ni Ni/Au, Au/Ni and Au/Au. (d) Scaled binding energies and separations for the cases 

shown In (a), (b), (c). 


Figure 7.— Variation of binding energy and force with separation for different bimetallic Junctions A-B (A.'tlp-B: substrate). 
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